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Abstract. Coulomb excitation is a standard method used to extract quadrupole 
excitation strengths of even-even nuclei. In typical analyses the reaction is assumed 
to be one-step. Coulomb only, and is treated within a semi-classical model. In this 
work, fuUy-quantal coupled-channel calculations are performed for three test cases in 
order to determine the importance of multi-step effects, nuclear contributions, feeding 
from other states and corrections to the semi-classical approximation. We study the 
excitation of ^°S, ^^Ni and '^^Kr on ^^'^Au at « 50 AMeV. We find that nuclear effects 
may contribute more than 10% and that feeding contributions can be larger than 
15%. These corrections do not alter significantly the published B{E2) values, however 
an additional theoretical error of up to 13% should be added to the experimental 
uncertainty if the semi-classical model is used. This theoretical error is reduced to less 
than 7% when performing a quantal coupled-channel analysis. 
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Nuclear collectivity of an even-even nucleus is closely related to its quadrupole 
electric reduced transition probability B{E2,0i —>■ 2^). This strength can be 
determined experimentally by measuring either the lifetime or the Coulomb excitation 
cross section of the 2^ excited state [H [2] . Originally, the Coulomb excitation technique 
(referred to as "Coulex" in the following) was used to measure properties of the target 
(e.g. OH]). The sub-Coulomb energies at which the reaction took place ensured a 
nuclear-free measurement. In the last decade, Coulex has been expanded to intermediate 
energies with the aim of studying unstable nuclei [2] . In this case the nucleus of interest 
is the beam particle and a heavy target is used to produce the virtual photons. The 
reaction takes place at high enough energy to inhibit multi-step effects and data is taken 
only at very forward angles, where one expects to be free from nuclear interference. This 
method has enabled accurate measurement of the B{E2) of a large variety of systems 

A systematic comparison between the intermediate-energy Coulex method and 
the lifetime method showed that there is consistency between the two techniques pj. 
Moreover, the accuracy of the B{E2) strengths extracted through Coulex is comparable 
to that from the lifetime measurements [1]. The work of Cook et al. [Ij focused on the 
experimental accuracy but did not consider the uncertainties due to the approximations 
in the Coulex theory used to connect cross sections and electric strengths. This is exactly 
the focus of the present study. 

One way to analyse these unstable beam Coulex experiments is using the semi- 
classical model of Alder and Winther [18]. It conveniently provides a linear relation 
between the Coulex cross section and the reduced transition probability. The 
approximations in the Alder and Winther theory [18] are three fold: 1) the straight- 
line semi-classical approximation; 2) the excitation is a one-step process; 3) it is 
purely Coulomb. The straight-line approximation is partially corrected within Alder 
and Winther [18]. The second point is not so straight-forward. Most of the nuclei 
studied through this technique exhibit large collectivity and thus have other excited 
states that are strongly coupled to either the ground state or the first excited 2"*" state. 
Even if it is generally assumed that the cross sections to these other excited states 
are small at intermediate beam energies, multi-step mechanisms and interferences can 
distort the desired result. In order to solidify the reliability of the intermediate-energy 
Coulex method, it is important to evaluate the uncertainties coming from the one- 
step approximation. Finally, nuclear contributions need to be consistently included 
in the calculations so that Coulomb-nuclear interference is correctly accounted for. 
The inclusion of nuclear effects in realistic quantum coupled-channel calculations may 
enhance multi-step effects. 

As mentioned above, intermediate-energy Coulex relies on restricting the scattering 
angles taken into account for integrating the cross sections to a range corresponding to 
impact parameters larger than the sum of the target and projectile radii. A detailed 
study of the sensitivity to the impact parameter cut was performed on the ^^Ar data 
[To] and results validate the procedure. In cases where low statistics forces the inclusion 
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of a wider angular range, nuclear effects have been estimated with quantum distorted 
wave calculations to be of the order of 6% [HI |15]. However, this value should not be 
taken as definitive, since, as we shall see, the nuclear contribution depends strongly on 
the particular analysis considered. 

Another problem that is considered when analyzing Coulex data is the possibility 
of feeding: the reaction process excites high lying states that could then decay to the 
2i state, producing an enhanced 2i Oi signal. In most of the studies, estimates 
of feeding predict it to be unimportant (e.g. [131 IE]) mostly because at intermediate 
energies the relative cross sections to higher spin states are small and larger excited states 
are hindered compared to the lower transitions. Nevertheless there have been cases 
where feeding needs to be carefully considered before a reliable strength is extracted 
[T2l [8] . In intermediate-energy Coulex experiments, the statistics is often low and the 
efficiency of the 7-ray detectors is limited such that a 7 peak for a feeding transition is 
rarely seen [8] . Feeding corrections are based on theoretical estimates [18] and subtracted 
from the 2^ cross section, before extracting the B{E2) strength. 

Intermediate-energy Coulex has been applied mostly to intermediate mass nuclei 
bound by a few MeV but as beam intensities improve, it will be applied to more 
exotic systems. The loosely bound nature of unstable nuclei has modified many of 
the traditional views of nuclear reactions. For example, when the exotic nucleus has an 
extended tail in its wave-functions, one finds nuclear contributions at impact parameters 
much larger than the sum of the target and projectile radii [19]. In addition, due to 
the proximity to the continuum, multi-step breakup effects need to be considered |20j . 
A systematic study of nuclear interference in the Coulomb dissociation of halo nuclei 
has shown large nuclear effects, even in the forward angular regions considered safe for 
Coulomb experiments [21]. A recent comprehensive study of the Coulex of ^^Be for 
extraction of the B{E1) between the two bound states validates the Coulex method 
across a wide range of beam energies, provided all these effects are taken into account in 
the theoretical model [22]. For the B{E2) of intermediate mass nuclei, it is important 
to solidify the theoretical methods used at present before these new dripline challenges 
can be faced. 

Table 1. Information on the intermediate-energy Coulex experiments considered here. 
For each case we give the beam laboratory energy, the maximum centre-of-mass angle 
for cross-section integration, the corresponding cross section for the 2^ state and the 
B{E2,0^ 2^) value extracted through Winthcr and Alder's theory [15] , 



Nucleus 


Energy 


nmax 


^1 


i?(i?2,0+^2+) 




(AMeV) 


(deg.) 


(mb) 


(e^fm^) 




35.7 


4.56 


39.6(3.8) 


350(33) 


58Ni 


72.4 


4.26 


175(36) 


707(145) 


78Kr 


57.4 


4.24 


1124(133) 


6244(738) 
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Table 2. Spin, parity and excitation energy for all the states included in the coupled- 
channel calculations. The O''" assignment for the 3.666 MeV state in ■^"S is based on 
a comparison of the experimental spectrum with the spectrum of the mirror nucleus 
and a shell-model calculation for the A = 30 isobars. 
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2t 
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2.775 
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4^ 
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9+ 
^3 
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In this work we perform fully-quantum coupled-channel calculations for three test 
cases that have been measured by intermediate-energy Coulex: ^°S, ^^Ni and ^^Kr. 
These three test cases span a variety of physical situations. The first, ^°S, corresponds 
to a very short lived isotope, two nucleons away from the proton dripline, with only a 
few excited states. The Coulex of ^^S was measured at 35.7 AMeV on ^^"^Au [9]. The 
second, ^^Ni is less exotic, contains very strong transitions to higher energy states and 
therefore has an important feeding correction. It has been measured several times before 
and we consider here the experiment at 72.4 AMeV on ^^''^Au [12]. The third, "^^Kr, has a 
very small 2^ excitation energy, and consequently a very large B{E2). Here we consider 
a recent measurement at 57.4 AMeV on a ^^^Au target [11]. Experimental details for 
these experiments are summarized in Tabled], where we include the B{E2) extracted in 
the corresponding studies using the first order semi-classical theory |l8j . 

We have investigated the spectra of these nuclei in detail and isolated the states 
that can affect the reaction mechanism. These are summarized in Table [2l For the two 
heavier cases, the spectra are well known. However, for ^°S, the spin and parity of the 
3.666 MeV state are undetermined. We have assumed they are O"*" by comparison with 
the level scheme of the ^°Si mirror nucleus [23] and a shell-model calculation for the 
T = 1 states of the A = 30 isobars [23]. As it can feed into the 2f state, it needs to be 
included in the calculations. 

States with unnatural parity (e.g. 1^ and 3^ states) were not included since they 
would decay to 2"^ and O"*" states by magnetic transitions which are not implemented in 
our coupled-channel calculations. In the energy range of interest, there is one l"*" state 
in ^''S, one 1"*" state in ^^Ni and none in ^^Kr. Generally, magnetic transitions in Coulex 
are much weaker than the electric ones [25] . 

A very large amount of transitions is possible between the excited states listed 
in Table [2l In Table [3] we list all the transitions taken into account in our coupled- 
channel calculations. We also provide halflives and branching ratios [23j from which we 
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Table 3. Transitions included in the calculations: spins and parities of the states 
involved in each transition, the corresponding halflives branching ratios /-,., the 

reduced transition matrix elements M{EX) and deformation lengths S = (3R. Unless 
otherwise noted, the spins, parities, halflives and branching ratios were taken from the 
NNDC database [53]. 





TTT 
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2| 


o| 




1.00 




18.71 " 
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3.29 


0.231 




2^ 


2^ 
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0.80 
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2.842 
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2| 


1 ps 


1.00 




9.31 


0.292 


c;q TV T- 

58^1 


2^ 


^ -1- 
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26.59 " 
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0^ 


970 fs ^ 


2x10"^ 
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0.718 




J J- 

4^ 


2^ 


970 fs ^ 


1.00 




71.56 


1.031 




2^ 


0^ 


0.38 ps 


4.3x10" 


_2 


1.40 


0.045 




2^ 


2^ 


0.38 ps 


0.96 




42.13 


0.607 




r-v-l- 

2^ 


4^ 


0.38 ps 


5.7x10" 
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36.75 


0.395 




2^ 


0^ 


52 fs 


0.40 




9.18 


0.296 




2^ 


2^ 


52 fs 


0.58 




56.58 


0.815 
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2^ 


A -X- 
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52 fs 
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49.62 


0.533 
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52 fs 
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72.18 


0.775 
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35 fs 


0.59 
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0.367 
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0.39 
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1.040 
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1.00 




79.02 


1.793 




0+ 


2| 


5.8 ps 


1.00 




41.89 " 


0.425 " 




4^ 


2| 


2.5 ps 


1.00 




125.68 


1.275 






o| 


3.7 ps 


0.39 




12.21 


0.277 




9+ 


2| 


3.7 ps 


0.61 




54.29 


0.551 



" using B{E2, 0^ ^ 2^) measured by Coulex P [HI [H] 
experimental lower value 

of the 4^^ 0^^ transition in "^^Ni 
from the 2X ^2 transition in ^^Ni 
assuming S(E2,0^ ^ 2|) = B{E2,A\ 2^) 



determined the B{EX). Reduced matrix elements M{EX) were evaluated directly from 
B{EXy. 

M{EX, h ^ h) = ^J{2h + l)B{EX,h^ h). 
For the Oi 2f transitions, we used B{E2) directly from the Coulex experiments 
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Table 4. Onc-stcp (DWBA) and coupled-channel (CC) calculations including only 
the ground state and the 2^ state: comparison between cross sections integrated over 
the experimental angular range for the full process with those from Coulomb only. All 
cross sections are in mb. 





DWBA(2^) 


CC(2^) 


30 S Coul+nucl 


44.7 


43.6 


=^°S Coulomb 


40.6 


40.6 


s^Ni Coul+nucl 


161.5 


161.1 


^^Ni Coulomb 


183.8 


182.7 


^^Kr Coul+nucl 


1052.6 


1030.7 


^^Kr Coulomb 


1056.3 


1036.0 



data under scrutiny [9l [121 E] • We have not included transitions with reduced matrix 
elements M{EX) two orders of magnitude lower than the main 0^ —> 2f transition. 
Since our coupled-channel calculations are restricted to electric couplings, transitions 
involving a mixture of Ml and E2 transitions (such as 2+ 2^ transitions) were 
assumed to be 100% E2. This is not always a good approximation. However, none 
of these transitions is involved in the direct excitation of the 2\ states, but in feeding 
couplings. Our aim was to study the effects of couplings, therefore we chose to keep the 
full strength of these transitions. In our calculations we assumed all coupling amplitudes 
have the same sign. We have checked this assumption and verified that uncertainties 
due to sign changes are neghgible. 

A specific remark is needed for ^^Kr, as the lifetime of the 02" state is unknown. 
This state, as well as the ^ and 23" states, lies at about twice the excitation energy 
of the 2\ state. This suggests that these O2 , 2^ and ^ states are quadrupole 2- 
phonon states, with the 2\ being the 1-phonon state. In a pure vibrational model, 
B{E2,02 2i) = B{E2,Ai 2f) [26j. This is what we assumed in order to obtain 
the transition matrix element for the O2" 2f transition in ^^Kr. 

As part of our aim is to check the importance of nuclear effects, nuclear couplings 
were also included in the calculations. We have assumed the matter deformation to be 
the same as the charge deformation for all three cases. As these systems are not halo- 
like, this approximation should be adequate. The corresponding deformation lengths 
are also included in Table O For a given transition, the nuclear and Coulomb couplings 
were considered in both directions, but reorientation couplings were not included, as 
would be the case in a vibrational model. Optical model parameters were taken from 
elastic scattering studies. For the ^^S case, we used the optical model parameters for 
^°Ar + 2°^Pb at 44 AMeV [27]. For ^^Ni and '^^Kr on ^^"^Au, we used the parameters from 
an elastic scattering study of ^^Kr + ^°^Pb at 43 AMeV [28]. A value of rc = 1.2 fm was 
used for the Coulomb radius parameter. Coupled-channel calculations were performed 
using the code FRESCO ^29]. 
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Table 5. Full coupled-channel calculations including all transitions specified in Table 
[31 /j„_^2+ ^^'^ branching ratio for the J,^ 2^ transition [5^. a feed are the cross 
sections for feeding into the 2^ state. The value for a feed in bold is the sum of the 
cross sections contributing to cr{2f) for each case. 









Full CC 




Nucleus 


State 


a 
(mb) 


I Ttt .9+ 


feed 

(mb) 




2^ 


43.6 




45.5 




2+ 


2.2 


0.80 


1.8 




0^ 


0.1 


1.00 


0.1 






155.6 




188.3 






8.6 


1.00 


8.6 


58Ni 




1.9 


0.96 


1.8 




2+ 


19.1 


0.58 


11.1 




2t 


29.0 


0.39 


11.3 






1013.9 




1041.9 




ot 


2.8 


1.00 


2.8 






8.0 


1.00 


8.0 




2t 


28.1 


0.61 


17.2 



For an accurate calculation of the Coulex cross section at forward angles one needs 
to be very careful with convergence. Partial waves up to Lmax = 3000 — 6000 and a radial 
integration up to Rmax = 300 fm were needed in order to get a converged integrated 
cross section within the experimental angular range. Checks of the sensitivity of the 25*" 
excitation cross section to the optical model parameters were performed and we estimate 
an uncertainty smaller than 5%. Our calculations show negligible sensitivity to the 
Coulomb radius rc- It is important to note that the calculations we have performed are 
non-relativistic. Relativistic kinematical effects have been studied within the context 
of the first-order semi-classical theory ^30]. We have estimated the effect of relativistic 
kinematics by repeating the calculations at a corrected beam energy and obtained for all 
three cases modifications in the cross sections smaller than 5%. Therefore, if we neglect 
uncertainties in the transition strengths used, the theoretical error on our cross sections 
should be smaller than 7%. This should be added to the experimental errors. 

Results including only the ground and 2^ states are presented in Table HI 
We present 1-step (DWBA) and coupled-channel (CC) calculations. Cross sections 
were integrated over the angular range of the particular experiment. We performed 
calculations including both nuclear and Coulomb excitations in the transition matrix 
elements and we compare them with the results obtained with only Coulomb. In all 
three cases, O^*" ^-T' 2^ multi-step effects are small. Nuclear effects are of the order of 
10%, 12% and 0.4% for ^os, ^^Ni and ^^Kr, respectively. These are still within the 
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For and ^^Ni, nuclear effects are not negligible and indicate that the maximum 
angle used in integrating the experimental cross section should be carefully chosen. The 
contribution of the nuclear part of the interaction to the cross section is indeed very 
sensitive to this maximum angle. For ^^Ni, our tests show that decreasing the maximum 
centre-of-mass angle by only 0.5 degree cuts the relative nuclear contribution by a factor 
2. 

One can also compare our Coulomb-only DWBA cross sections to those predicted by 
the semi-classical model (see Tabled]). From this comparison we find that the straight- 
line trajectory approximation alone introduces an error in the cross section of 6% at the 
most. 

Full coupled-channel results including the states in Table [2] and transitions from 
Table [3] are presented in Table [5l The cross sections to individual states (a) are 
multiplied by the branching ratios to the 2^ state in order to get the feeding 
contributions. The sum of these and the Coulex 2f cross section gives the full 
cross section (in bold) to be compared to the experiment. Feeding contributions are 
important in ^^Ni (^ 17%), as we already knew, but also in (~ 5%). For ^^Ni, 
the feeding correction estimated in [12] of 25 mb is significantly smaller than our 
prediction (33 mb), mainly because of the contribution from the 4^ state which was 
omitted in the experimental estimation [12]. For our theoretical cross section taking 
B{E2, of 2i) = 350 e^fm^ is above the upper limit of the experimental range. A 
reduction of the B[E2) value to 303 e^fm^ is necessary to get a theoretical cross section 
in agreement with the experimental result. Considering the 7% theoretical error, this 
value is still compatible with the value extracted in [9j. For ^^Ni and ''^Kr, the expected 
cross section is within ^ 7% of the mean experimental value, i.e. within the experimental 
range. 

In conclusion, full coupled-channel calculations for ^"S, ^^Ni and ^^Kr Coulex 
confirm the agreement between B{E2) extracted through the lifetime and the Coulex 
methods [1]. Our study shows that theoretical contributions to the errors in the cross 
sections can be ~ 13% if the first order semi-classical theory of Alder and Winther is 
used [18] but should be less than 7% if a full quantum coupled-channel calculation is 
performed and feeding is consistently taken into account. 

We thank Alexandra Gade for useful discussions and comments on earlier versions 
of the manuscript. This work was supported by NSCL, Michigan State University, and 
the National Science Foundation through grant PHY-0555893. 
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